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We present the calculation of the density matrix response function of the two-dimensional electron 
gas with Rashba spin-orbit interaction, which is applicable in a wide range of parameters covering 
the diffusive and non-diffusive, the dirty and the clean limits. A description of the crossover between 
the different regimes is thus provided as well. On the basis of the derived microscopic expressions we 
study the propagating charge and spin-polarization modes in the clean, non-diffusive regime, which 
is accessible in the modern experiments. 
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I. INTRODUCTION 

One of the important directions of study in the field of 
spintronicsi is a possibility to manipulate spin-polarized 
electrons by means of electric fields. In two-dimensional 
heterostructures or quantum wells, this can be achieved^ 
due to Rashba spin-orbit (SO) interaction'^ which arises 
from the structure inversion asymmetry and couples lin- 
early the electron's momentum to its spin. 

The standard approaches to the study of the charge 
and spin dynamics in impure two-dimensional electron 
gas (2DEG) with Rashba SO coupling is based on either 
Kubo formalism."* or quantum kinetic equation^ (QKE). 
In the diffusive limit vpq,u} ^ , where r is an elas- 
tic mean free time, the usually made approximation is 
a small-g expansion of the so-called diffusion propaga- 
tor (or, equivalently, the angle averaging of the QKE). 
As a result, one ends up with the diffusion equation,i^ 
which along with the Dyakonov-Perel spin relaxation 
mechanism,^ accounts for the spin precession as well as 
for the dynamic spin-charge mixing. 

Such description applies at the large distances and time 
scales. The study of the charge- and spin-density evolu- 
tion on short [vfI,^ S> t^^) or intermediate [vfI,^ ^ 
T~^) scales requires either a more refine consideration 
of the density matrix response function or a solution of 
the QKE in the most general case. The recent progress 
in experimental observation of the non-diffusive spin 
dynamics^ in the transient spin grating experimentsSiS 
makes desirable the further elaboration of theoretical 
tools applicable beyond the diffusive limit. 

While comparing the above mentioned approaches, 
we would like to remark that the QKE is the differen- 
tial equation for the local, non-stationary distribution 
functioniS (more precisely, for the Keldysh component 
of a Green's function), which in the presence of the SO 
coupling has been recently considered in the spatially 
nonuniform situations in Refs. [Ill and [H. In the quasi- 
equilibrium regime the QKE approach can be comple- 
mented by the alternative, but equivalent, consideration 
of the charge and spin densities in the Fourier represen- 
tation. This seems to be especially useful in the context 



of the boundary problem, which appears to be a subtle 
issue of the QKE approach in the clean, non-diffusive 
regime^^. 

In this paper we present the compact analytical ex- 
pressions for the density matrix response function of the 
impure 2DEG with Rashba SO coupling. They cover 
quite a wide range of parameters, extending from the 
scale of the diffusive transport to the arbitrary spatial 
and temporal modulations of external fields as well as to 
the arbitrary values of the cleanness parameter aj^rkp, 
an being the Rashba constant. On this ground the de- 
scription of crossover from the diffusive to non-diffusive 
(ballistic) regimes becomes quite efficient. In particular, 
we apply the obtained expressions to the study of prop- 
agating charge and spin-polarization modes, and discuss 
the (/-dependence of their dispersions and lifetimes. 



II. CHARGE AND SPIN DENSITY RESPONSE 
FUNCTIONS 

The Hamiltonian of the 2DEG with Rashba SO 
coupling^ in the absence of impurities reads 



H{k) = — + aii{(j''ky ~ ayk^), 
2m 



(1) 



where m is an effective mass, and fi, = 1 in the units 
used. The spectrum of ^ is given by — + fj,aii\k\, 
where /i = ± is a subband index. The projectors onto 
the corresponding eigenstates are found to be 
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Using ([2]), one can perform a spectral decomposition of 
the Hamiltonian ([T]): 



(3) 



A weak disorder is modelled by random isotropic point- 
like scatterers with Gaussian configuration distribution, 
and its strength is parametrically given by (mr)"*. 
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For the future reference we introduce the following no- 
tations for the Fermi momentum kp — y^2mEp + kj^, 
the Fermi velocity vf = and the Fermi energy Ep; 
as well as for the Rashba momentum splitting kji — man 
and the density of states per spin component — j^- 
When appropriate, we will also use the dimensionless 
units y = z = ^> = M = i§r{u^ + ^) ^ ^(i. 

Like in Ref. [l^ we can define the following regimes 
depending on the relative strength of SO coupling and 
disorder potential: dirty y <^T, clean <T <^ y, and 
super-clean F < y^. The last regime is hardly accessible 
experimentally by now, and it will not be discussed in this 
paper. However, it is important to note that the clean 
regime as it is introduced here implies the limitations on 
F from below as well as from above. The necessity to dis- 
tinguish between the clean and the super-clean regimes 
will be discussed in the following. 

A linear response of a system p to an external pertur- 
bation V is determined by the density response func- 
tion X = x(q, w). In case of spin-resolved spectrum, 
X = (x"'') is a 4 X 4-matrix {a, (3 = 0,x,y,z; the in- 
dex "0" corresponds to the charge component), which 
includes the polarization operator (a,/3 = 0), the mag- 
netization (a,/3 = x,y,z), and the spin-charge mixing 
(a = 0, (3 = x,y,z; or vice versa). Quite generally, x 
is evaluated by means of the disorder-averaged Green's 
function technique^'*. In the self-consistent Born approx- 
imation one finds at zero temperature^ X = + X^^j 
where 

x\ci,uo)^2ivT / d£/(q,c^;e)[l~/(q,c^; £)]-!, (4) 

^^.'C-'^+q) (5) 

{Ep + s-^~e^){Ep+e + ^-e'i^^^y' 

and .F;^^,(k,k + q) = TrKP^(k)a"P^,(k + q)], (7° = 1, 

u) = uj -\- - , as well as x^^''^^ = 2v5°'^ . A careful analysis 
shows that 

X\ci,uj) « 2^l.c^T/(q,c^)[l - /(q,c^)]-\ (6) 

where /(q, w) = /(q, a';£ = 0), remains a good approxi- 
mation in the dirty (y ^ F), intermediate {y T), and 
clean (y^ < F <C y) regimes. It breaks down only in 
the super-clean regime, y^ > F, and in this case the in- 
tegration over £ in Eq. ^ has to be performed exactly. 
On the other hand, in the super-clean limit the disorder 
potential can be neglected at all (F ^ 0). Then, follow- 
ing the approach of Ref. [3, one can find exact analytic 
expressions for x"'^ at arbitrary q and lo in this regime. 

The vertex corrections to x are accounted in (|4]) and 
^ in the diffusion propagator, or diffuson, 'D{(\,uj) = 
t[\ — /(q, cij)]"^. In the diffusive regime vpq,bj <C t^^, 
Eq. ^ can be expanded in a series of small q and u, 



which makes straightforward the evaluation of /(q, lS) 
and I?(q, a;)'*. The Fourier transform of 2?~^ gives rise 
to the diffusion equation^i^ describing the evolution of 
charge and spin densities at large distances and time 
scales. In turn, a consideration of the short-range dynam- 
ics requires a more careful evaluation of the momentum 
integral in ([5]), especially of its angular part. 

Before proceeding further in evaluation of x, let us es- 
tablish the relation between the discussed diagrammatic 
approach and the approach based on the quantum kinetic 
equation^. Both of them can be deduced from the Dyson 
equation- . We refer to the Appendix \K\ where the corre- 
sponding alternative derivation of the response functions 
dU,® is presented. Note that the quasiclassical condi- 
tions q <^ kp and co <C Ep have not been applied yet for 
the evaluation of x- This will be done in the next section. 

In turn, the derivation of the QKE from the Dyson 
equation is based on the procedure which is specifically 
designed for the quasiclassical limit^*^. The QKE rep- 
resents the differential equation for the nonequilibrium 
distribution function. In the quasi-equilibrium (linear re- 
sponse) regime the QKE approach becomes fully equiv- 
alent to the quasiclassical linear response theory, since 
both are derived from the Dyson equation under the same 
assumptions. 

The solution of the QKE for the system in question 
also allows us to access the charge and spin dynamics in 
the clean and dirty, diffusive and non-diffusive cases. In 
the spatially nonuniform and non-diffusive situations the 
QKE in presence of the SO coupling has been previously 
considered in Refs. [ll| and [11. However, it is important 
to note that the time and space representation of the lo- 
cal distribution function requires a proper choice of the 
boundary conditions, which in presence of the strong SO 
coupling and weak disorder becomes a nontrivial prob- 
lem. For example, in Ref. [l^ the boundary conditions 
imposed on the distribution function in the clean, non- 
diffusive case have been motivated by the phenomeno- 
logical arguments borrowed from the scattering matrix 
theory, on the analogy with Ref. [H. We also refer to 
Ref. llTI . where similar issues have been analyzed in the 
diffusive regime. 

The problem of the boundary conditions can be post- 
poned (and, hopefully, circumvented), if the response 
functions are evaluated in the frequency and momentum 
representation. For example, applying the inverse Fourier 
transformation to a finite system of the size L ^ kp^ , we 
would obtain the spatial dependence of the observable 
quantities p"(r,t). The boundary conditions are thus ef- 
fectively accounted in terms of y"(r,<). 



III. QUASICLASSICAL APPROXIMATION 

The quasiclassical approximation of x and / [Eqs. (|4]) 
and ^] is achieved due to the substitution of the integral 
over rf|k| in Eq. (I5|) by a residue value. To this end, 
one can close a contour of integration in either upper 
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or lower half-planes. However, a better approximation 
would be to take the mean value of the both results, 
because it supports the symmetry /^(q, = /(— q, —lo) 
of the exact expression ([S]) (at £ = 0). Thus, we obtain 
the sum of the residia 



ra/3 



ikX 2 



87rr(fe + iikn) 



^:f'(k,k + q) 



Aw 



'k+q 



(7) 

calculated in the poles fc^ ^ = V^f ~ X'm*Cj — ^kji. The 
new index A labels the two different closures of the in- 
tegration contour. Note that we also assume <C kp^ 
which is a physically justified condition. 

In the basis q||ea:, the integration angle 4> is mea- 
sured from the direction of q, and all non- vanishing com- 
ponents (sec the Appendix |B]) are distributed into the 
two blocks: (I^o jyyjoy ^ jyO) ^^^.^ ^^le even block 
(a,/3 = 0,y), while (/^^^/^^^/^^ ^ _/^^) form the odd 
block (a,/3 = x,z). The parameter s reflects the sym- 
metry of JF"^, (k, k + q) = sJF",'^^(— k — q, — k) under the 
three successive operations: k ^ k — q, k — > — k, and 
/i ^ /i'. It appears to be that s equals 1 for all compo- 
nents, except for s = — 1 for 7°^ = 

The central point of our calculation is an exact inte- 
gration over the angle variable in Eq. ([7|). The technical 
details of this evaluation are described in the Appendix 
IB] Here we present the final result for the non-zero com- 
ponents of /"^: 
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r(e) 



JZX JXZ 



jyy 



= I. 



where 



Tie) _ Sr 'r(-M)" ( C 



j(o) _ »r(-M)" ( [T^ 



(8) 
(9) 

(10) 
(11) 

, ,u ) and 

n = 0,1,2. Note that li^^ and In"' are mapped onto 
each other, J^'^'' 



In these expressions = — , C — Cn 



when C C/j a-nd z ^ are 
under the formal 



exchanged. Moreover, In'°^ —In'°^ 
replacements C, ^ z and z^ ^ C/j- 

The expressions (|8|)-([TT|) represent the main result of 
this paper. Their applicability range is determined by 
the conditions y,r, <C 1 and <C F which follow 
from the assumptions made in the course of the calcu- 
lation. As for the rest, they are equally well applicable 
in the diffusive {z^w <ti F) and non-diffusive {z,w S> F) 
regimes, in the dirty (y <C F) and clean {y 3> F) limits. 

It is useful to consider different limiting cases of 
Eqs. ([8|)- pT|) . in order to prove their consistency with 
the well-established results. 

In the limit y — * 0, the off-diagonal terms vanish, 7°^ — 
JZX _ -iffl^ile the diagonal terms yield the familiar 



expression in the absence of SO coupling, 



^y=o - 



iV 



(12) 



In the homogeneous limit z 



0, = C, / 

' in ' 



and /^^ = = i(/°" + /^^). Setting w = in 
these expressions, we recover the out-of-plane — 7(1 — 



and the in-plane 



^ > " r[i+(y/rp\ 
^(1 — = 2^ spin relaxation rates of the Dyakonov- 

Pcrcl mechanism^. In turn, the small-z behavior of the 
off-diagonal components reads 



iy^Tz 



2w^(w^ — 



JZX 



yTwz 



7,2 



y 



2^2 ■ 



(13) 



At ui = 0, we recover the coefficients \imz^o[ipy / (2fci?z)] 
and \im2^a[il'^^ /{2kpz)] of the diffusion equation de- 
scribing the spin-charge mixing and the spin precession, 
respectively, which were previously derived in Rcfs. 4.^. 

We note that the exact integration over the angle vari- 
able yields the expressions (jTU)) and pT|) . which contain 
the square roots of the complex-valued expressions with 
the branch-cuts along the lines w = z and w = \z ±y\. 
In the next section we will show that these branch-cuts 
determine dispersions of propagating modes in the clean, 
non-diffusive limit. Meanwhile, the perturbative expan- 
sion in z, which is only validated in the diffusive limit, 
destroys the square-root structure of the expressions (fTUl) 
and (jlip . thus eliminating the branch-cuts from them. 
Therefore, one can hardly expect an occurrence of dis- 
persive propagating modes in the diffusive regime. 



IV. STRUCTURE FACTORS AND 
PROPAGATING MODES 

Diffusive modes and their lifetimes^ still remain in fo- 
cus of an active study. They are naturally reproduced 
within the present approach in the dirty limit y ^ F at 
z, w <C F. However, there exists another interesting ques- 
tion about a probable occurrence of propagating modes 
in the clean regime F ^ y. For example, in Ref. [l^ 
it was conjectured on the basis of quite general physical 
arguments that "in the limit of strong spin-orbit cou- 
pling, there is a regime where a propagating, coupled 
spin-charge mode is possible." Now, the microscopically 
derived expressions (|5|)- (fTT|) allow us to tackle this is- 
sue as well as the issue of propagation of coupled s^-s^- 
polarization modes determined by the odd block of x- 

A convenient tool for studying the dispersions and 
lifetimes of either propagating or diffusive modes in all 
regimes is provided by the structure factor 5(q, w) — 
Tr Imx(q, w), which is a directly measurable functioni^. 
In case of Rashba SO coupling we can define the two 
block-resolved structure factors, S*^*^' = Trelm^ and 
= Trolmx, such that S = S^"^ + S^°\ The peaks' 
positions of 5'('^'°) in the plane {q, lj) correspond to the 
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FIG. 1: The structure factors S^''\z,w)/u (upper panel) and 



S^°'>{z,w)/u (lower panel) at y = 10"^ and T = 0.5 x 10" 



dispersions of propagating (at w = - diffusive) modes, 
while the peaks' widths determine the inverse hfetimes 
of the corresponding modes. 

In the upper and lower panels of the Fig. [T] we plot 
S'^^^z,w)/iy and S'-°^z,w)/v, respectively. The SO and 
the disorder parameters are chosen to be y = 10~^ and 
F = 0.5 X 10"''. By the order of magnitude these val- 
ues correspond to those of the high-mobility sample from 
Ref.'?' which was characterized by /x « 1.5 x 10^ cm^/V-s. 
Note that y ^ F ^ = 10~^, i.e. the chosen parame- 
ters belong to the clean regime, but they are still far away 
from the super-clean regime. The values of the grating 
wave vector— q = 0.44 — 5.3 x lO'' cm~^ correspond to 
z = 10^^ — 10~^, which means that the non-diffusive 
regime z > F is accessible. However, a direct comparison 
with the experiment is obstructed by an admixture of 
the Dresselhaus contribution to the SO coupling, which 
is not taken into account in our consideration. 




FIG. 



Upper panel: the peaks' heights (weights) 



5'<°'(z,u>£°'(2)) (dashed and dotted) and S^''\z 



(solid) normalized by u. Lower panel: the widths ■y± \z) and 
7+ '(2) of the corresponding modes normalized by F. 



In the Fig. [T] one can clearly see the dispersions of the 
four propagating modes: w^jf'* ~ z and w^^-* « y in the 



upper panel, and w^f'' « z -I- y and w 
the lower panel. We note that the fact of existence of 

(e) 

the gapless mode w\_ « z has been already pointed out 
in Ref. To' as the consequence of the charge conservation 
law. 

In order to understand physical meaning of the other 
modes, let us consider S^'''°\z, w) at small and large z. 



' y 111 
z — y\ in 



In the homogeneous limit z 



0, w 



(e) 



,(°) 



y, 



and we can approximate the structure factors by the 
Lorentzians, 



S'(^'°)(0,w) 



2iyy 



(e.o)i^\ (e,o) 



(0) 



(14) 



with the weights a^_^\o) = ^F and a^^\Q) = ^F, and 

the widths 7i"^(0) = |F and 7i°^(0) = ^F. 

With increasing z the mode w'f'^ broadens very fast, 
simultaneously losing its weight. At the value z « 0.5y 
it is no longer resolved, but after this point its weight 
is augmented again. Thus, the saddle near the point 
{z,w) ~ (0.5?;, y) is formed in the landscape of S'^'^\ 
At z ~ y both and w^'' merge together, and 

for large z ^ y the mode w'^'' becomes degenerate. 
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in the sense that the dispersions of propagating charge 
and s^-polarization modes overlap. It also means that 
= Im(x™ + X^^) ~ 2Imx°" « 2S'o(z,w), where 
So{z,w) = ^''^^™-\J' w-z+iT ^ large-z structure factor 
in the absence of SO coupling. 

We note that even in the clean limit the value of /"^ 
remains negligibly small (|/"''| ^ |/°"|,|/^^|) at any z 
and w. By this reason, the spin-charge mixing can hardly 
be detected in the study of the propagating modes. 

The modes w''^^ and w''"^ behave more peculiarly at 
finite z, especially the latter one. Both of them are cou- 
pled s^-s^-polarization modes originating in the point 
(Zjw) « (0,?/). However, there is a qualitative differ- 
ence between them: the mode w^'' grows monotonously 

with z, while the mode idl' initially goes down towards 
{z,w) K, (y, 0), where it turns up, and after that it par- 
allels w''^\ The asymptotic behavior of S^""^ at large 
z ~ w > ?/ is S''°\z, w) « Sq{z, w - y) + S'o(z, w + y), 

i.e. Wj?' remain well-resolved and split-off of vJ"^^ . Note 
that the SO parameter y determines only the amount of 
this splitting, but not the shape of the individual peaks. 

The most interesting aspect seen in the lower panel of 
the Fig. [1] is the form of the structure factor 5*^°^ in the 
vicinity of {z,w) « (y, 0): it locally vanishes, which cor- 
responds to the local suppression of the w'^°^ mode. This 
is also clearly manifested in the upper panel of the Fig. [21 
where the absolute values (weights) of the structure fac- 
tors 5'('=)(z,z) and S'^°\z,\z ± y\) are plotted. Using 
the Lorentzian approximation, we estimate the widths 
7^^'^^(z) of the corresponding peaks, and present them in 
the lower panel of the Fig. [21 One can notice that near 

z ~ y, ')^°\z) (dotted curve) shows up specific features 
which are, however, somewhat artificially created: since 
the corresponding peak locally vanishes, the Lorentzian 
approximation also becomes locally inadequate. 

Thus, we have shown that the propagating spin-density 
modes may occur in the non-diffusive regime, provided 
the SO coupling parameter exceeds the disorder broad- 
ening. In order to demonstrate that the latter condition 
is important we also plot in the Fig. [3l the both structure 
factors for the values y = V = 10"'^. One can observe 
that now the peaks in the landscapes of 5*'-'^'°^ are no 
longer sharply pronounced and well-resolved as compared 
to those shown in the Fig. [1] The trend for decreasing 
ratio y/T now becomes quite transparent: the resolu- 
tion between the dispersions of the charge- and the spin- 
density modes is disappearing, and it completely vanishes 
in the dirty limit y/T<^ 1. 



V. SPIN-DENSITY RESPONSE TO 
INHOMOGENEOUS LONGITUDINAL 
ELECTRIC FIELD 

It is natural to expect that the landscapes of the com- 
ponents x''^? X^^: and x^" would manifest qualitatively 




FIG. 3: The structure factors S^'^\z,w)/u (upper panel) and 
S'^°Hz,w)/i' (lower panel) at y = T = 10"^. 

the same features as 5^*^^ while the landscapes of x^^, 
X^^ , and x^^ would have much in common with S^°\ 

In this Section we briefiy discuss the finite-g behav- 
ior of the off-diagonal spin-charge component x^°- It is 
linked with s*'-response MJf to inhomogeneous longitudi- 
nal electric field Eq^^ = —igV^^ applied in the x-direction: 

^"^^ " " ~4fc7rz[(i - /oo)(i - lyy) - {1^^)^] ' 

(15) 

In the homogeneous limit we obtain the expressio n -'^^i^^ 

limAf^ = — — - — — ^3 — ^-ir- — ^. (16) 

9^0 ^ 2kF 2w{w^ - y^) - iTy^ ^ ' 

In the clean regime y Y Eq. ([16]) exhibits the pro- 
nounced resonance peak at « y. Taking the stationary 
limit w ^ 0, one also recovers from the result of 
Ref. HO. 

The inhomogeneous response function (|15p is expected 
to possess interesting features in the clean limit as well. 
In the Fig. [H we plot the normalized as a function 
of the frequency at different values of the momentum 
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FIG. 4: The real (solid) and the imaginary (dashed) parts of 

M« = {4kF/eu)My{y X {0.01 + 0.15j}, w) + 18j [j = 0, . . . 4] 
at y = 10"^ and T = 0.5 x 10"^ 

transfer. Solid and dashed curves correspond to the real 
and the imaginary parts of respectively. One can 
observe that with increasing z the homogeneous peak at 
w = y is getting suppressed, while there develops the 
other peak which has the opposite sign and is centered 
at w = z. 



VI. CLEAN VS SUPER-CLEAN REGIMES 

Let us now present a posteriori justification for distin- 
guishing between the clean < T < y and the super- 
clean r < regimes. 

As one can see from the expressions pUj) and dTTI) . 
in the absence of the disorder broadening (F 0) the 
response functions x would diverge at w'^'^ = z and 
w'^ — y± z, i.e. at the same positions where the disper- 
sions of the propagating modes in the clean regime have 
been observed in Sec. IIVI In the limit F ^ these lines 
can be associated with the intrasubband and the inter- 
subband boundaries of the Landau particle-hole excita- 
tion region which is defined by the condition Imx"" ^ 0. 
The fact that the boundaries in question have the linear 
form results from the quasiclassical approximation. 

On the other hand, in Ref. [l^ it has been established 
that the intersubband Landau region in the presence of 
the Rashba SO coupling is bounded by the parabolas 
w± = (?/ ± z) ± (y ± z)^, and the imaginary part of the 
polarization operator does not diverge near them, but 
rather smoothly decays to zero. We note that the ap- 
proach of Ref. Il5| did not involve any sort of quasiclassical 
approximation, and was based on the exact evaluation of 
defined in Eq. (IXTtI) . 

Thus, one can conclude that the quasiclassical approx- 
imation in the super-clean regime breaks down not only 



aX z,w ^ 1, but also near the boundaries of the inter- 
subband Landau region. In fact, one does need some 
amount of disorder broadening F to regularize x near 
the lines — y ± z, and to thus justify the validity of 
the quasiclassical approximation in their neighborhood. 
The quantitative criterion for this value of F is 

F>max|Aw±|, (17) 

where Aw± = w± — w'^ = ±(j/ ± z)^ is the mismatch 
between the actual and the quasiclassical Landau edges. 

Since we are mostly concerned with the situations 
z ~ y, we obtain the condition F > which defines 
the clean regime. In the super-clean regime T < the 
quasiclassical approximation is no longer applicable, and 
therefore the study of the response functions (jA17p re- 
quires a separate consideration similar to that of Ref. [l^. 



VII. SUMMARY 

We have presented the compact analytical expressions 
for the density matrix response function of the 2DEG 
with Rashba SO coupling calculated at arbitrary ratio 
of vpci, ajikp, and r~^. They describe the linear 
response of charge and spin densities to the spatially and 
temporally modulated fields in any regime of interest, 
namely: clean and dirty, diffusive and non-diffusive. The 
calculated expressions (|10l) and (fTT|) . which define the 
density-density response to the nonuniform fields, as well 
determine the vertex corrections to the density-current 
and current-current response functions, and should prove 
useful in evaluation of the latter in the nonuniform case. 

The structure factors containing the information about 
the quasiparticle modes have been considered. We have 
particularly discussed the dispersions and lifetimes of the 
propagating charge and spin-polarization modes occur- 
ring in the clean, non-diffusive regime, which is experi- 
mentally achievable by now. The spin-density response 
to the inhomogeneous longitudinal electric field has been 
studied in this regime as well. 

The useful discussions with Matthias Eschrig, Vladimir 
Gritsev, and Alexander Shnirman are gratefully acknowl- 
edged. The work was supported by the Deutsche For- 
schungsgemeinschaft (DFG). 



APPENDIX A: DYSON EQUATION 

A very general approach to the study of the quasi- 
particle dynamics is based on the Dyson equationi^ 

(Go-i-I])*G-l, (Al) 

where G^^ = idt - H{k) - V{r,t), and V{r,t) = 
J2p cr^V^{T, t) is a time-dependent, spatially nonuni- 
form and spin-dependent external potential. The Green's 
function G and the self-energy E are the matrices in the 
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Keldysh space 



G 



Gr Gk 
Ga 



Gr Gk 
Ga 



(A2) 



which include the retarded (R), advanced (K), and the 
Keldysh (K) components. Eq. (|A1|) can be equivalently 
rewritten as the set of equations 



(Gq ^ — '^R.a) * Gr^a ~ 1, 



(Gg ^ — Si?) * Gk — * Ga 



0. 



(A3) 
(A4) 



They should be also complemented by the certain choice 
of the self-energy S. In the self-consistent Born approx- 
imation the latter is determined by the self-consistency 
relation 



S(ri,ti;r'i,i'i) 



1 



G(ri,ii;ri,i'i)5(ri-r'i). (A5) 



There are different ways to treat the set of equations 
(|A3P - (|A5|) . One possibility is to apply the gradient ex- 
pansion, and to thus derive the differential equation for 
the quasiclassical Keldysh function Gk^. Since Gk is 
related to the distribution function, the solution of this 
equation determines the evolution of the charge and the 
spin densitiesi^. 

Another possibility is to treat the equations (|A3[) - (|A5|) 
perturbatively in V , and to derive the charge and the 
spin response functions. This approach is equivalent to 
the Kubo linear response theory, and in general it does 
not require the quasiclassical condition to be fulfilled. 

Let us expand G and S in a formal series of V: 

G«GW+G(i), S«S(o)+s(i). (A6) 

One can find that the equilibrium functions G*^"^ and S'"^ 
read 



G 



(0) 
R.A 



E - (k) ± 



2r 



^(0) _ * 



Gf = h{E){Gf~Gfl 
4°) ^ h{E){^f^^fl 



(A7) 

(A8) 

(A9) 
(AlO) 



where h{E) = 1-2 f{E) = tanh f{E) is the Fermi 

distribution function, and T is the temperature. 
The deviation from the equilibrium is described by 



G 



(1) 

K 



g'-j^Kv^g 



(0) 
K 



G%Kv* G 



(0) 



K 



(0) 



(All) 



The last term in this expression accounts the vertex cor- 
rections which are the same as obtained by the sum- 
mation of the ladder diagrams in the diagrammatic ap- 
proach. 



Let us now introduce 5^ — ^Tr (cr"G^^ ) and a% 



iTr ^cr"S^''^ , and rewrite (jAlip in the Fourier represen- 
tation. We obtain 



a cid I a, II , a, vertex 

9k = 9 k + 9 k + 9 k 



(A12) 



where 



g'^/ = Y: [h{E + DffS^ - KE ^)g-^] V^, 
9-/' = Y: [h{E - ^)gl'R ~ h{E + ^)g-A] 



a, vertex \ a. 3 

9k = 2^ 9ra'^ 



The terms which appear in the above relations read 



Q/3 

9r{a),r{A) 2 



= i^TrKF,,(k+f)^^P,(k-§)] 



.(A13) 



The projective representation of (|A13[) is due to the spec- 
tral decomposition of G^-*^ similar to ([3]) . 

In order to find ctJ^ let us now integrate over k. 
Introducing the definitions 



^FHA)MA) 



1 



Q/3 



^2t^Y9r(a),r(A) 



(A14) 



and taking into account the self-consistency relation 
(lASp . we obtain 

UK = h{E + u/2)[lAR-lRA]V 

+ hiE - u;/2)IrrV - h{E + lu/2)IaaV (A15) 

+ [h{E + uj/2) - h{E - oj/2)]IraV + IraOk, 

where the upper (spin and charge) indices are omitted for 
brevity. From this expression we can find (JKi which is 
proportional to V , and consequently the linear response 
X of the system: 



■ j 'Sl-K =Yxf.yl- (A16) 



One can show that the term in the first line of (|A15|) is 
negligible, the second-line terms approximately give the 
constant term i ^-nd the third line determines after 
the identifications e — E — Ep and Ira = I ■ Thus, at 
zero temperature the expression Q is recovered. 

In the super-clean limit t ^ <x> the self-energy and the 
vertex corrections disappear, and we obtain from (jA12[) 
the definition of the response functions 



A-quj 



E 



TrKP^(k)a"P^,(k + q)] 



w + iO + ej:-£C;q' 



(A17) 
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APPENDIX B: DETAILS OF THE ANGLE 
INTEGRATION 

Let us explicitly list the components of JF^^, (k, k + q) 
introduced in ([5|). Along with the diagonal terms 

= l[l±cos(0k+q-0k)], (Bl) 

(B2) 



we 



= i[l±cos(0k+, + 0k)], 
have the ofF-diagonal spin-charge 

Kl' = -^^M' = \ (/^ sin '^k + m' sin 0k+q 



mm' mm' 



i (/i cos 0k + /i' COS ^k+q) 



(B3) 
(B4) 



Kl' = -Kl' = ^MM'sin((/)k+q - 0k), (B5) 
and spin-spin terms 



mm' 



■^mm' ^ 



•^mm' " ^(A^sin(/)k - //'sin0k+q), (B6) 
■^MM' = -2(A^cos0k - M'cos0k+q), (B7) 



MM' / 



MM' 



One can 
relations 



= sin((/)k+q + 0k)- 

further elaborate these expressions 

, q + k cos 



using 



(B8) 
the 



COs(0k+q 

sin(0k+q 



|k + q| 
fc sin 



(B9) 
(BIO) 



where = 0k — 0q- 



Without the loss of generality, we can choose the basis 
defined by qlle^., or 0q = 0. This considerably simplifies 
our calculation due to the isotropy of the spectrum in 
the momentum space. In particular, one can observe that 
in this basis the components /°^, /^^, and /^^ would 
identically vanish, since the corresponding integrals in ([7|) 
have the form 



2ir 



sin F(cos 0) = 0. 



(Bll) 



Let us make the following observation which allows 
us to perform an exact angle integration in ([7]) for the 
remnant components. Changing the dummy index /i' — > 
/iX, where x = and replacing — > we find 

after the latter summation that odd powers of |k + q| 
disappear, and the integrand becomes a rational function 
of COS0 (cf. Ref. IT5I ). It only remains to expand it into 
simple fractions, and to apply the standard integral 



1 

2^ 



27r 



1 



l-h5cos0 V]~~62' ^^^"^^ 
where h is the arbitrary complex number. Thus, after 
tedious, but straightforward, calculation one can recover 

Eqs. (Uni) and (mi- 

We also note that the expression for the components of 
/ in the arbitrary basis can be obtained by the orthogonal 
transformation / = Oq/O"^, where 



Oq = 



/I 

cos(; 

sin (/ 
V 





— sint 
cos 







1/ 



(B13) 



and 0q is the angle of the vector q in the new basis. 
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